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C/2 , 

JP i Abstract The central problem of this chapter is temporal coherence of a three- 

dimensional spatially homogeneous Bose-condensed gas, initially prepared at finite 
temperature and then evolving as an isolated interacting system. A first theoretical 
j^ I tool is a number-conserving Bogoliubov approach that allows to describe the system 

^ . as a weakly interacting gas of quasi-particles. This approach naturally introduces the 

"sJ I phase operator of the condensate; a central actor since loss of temporal coherence is 

governed by the spreading of the condensate phase-change. A second tool is the set 
of kinetic equations describing the Beliaev-Landau processes for the quasi-particles. 
We find that in general the variance of the condensate phase-change at long times 
i-^ I t is the sum of a ballistic term °^ t^ and a diffusive term °^ t with temperature and 

^ . interaction dependent coefficients. In the thermodynamic limit, the diffusion coeffi- 

O ' cient scales as the inverse of the system volume. The coefficient of f ^ scales as the 

inverse volume squared times the variance of the energy of the system in the initial 
state and can also be obtained by a quantum ergodic theory (the so-called eigenstate 
thermalisation hypothesis). 



1 Description of the problem 



m 

f^ . We consider a single-spin state Bose gas prepared in equilibrium. To extract the 

CN ' relevant physics, we avoid the complication of harmonic trapping present in real 

experiments H] |2l [3] and we consider a spatially homogeneous system in a par- 
allelepipedic quantization volume V with periodic boundary conditions. In all the 
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chapter except subsection |4.3| the total particle number is fixed and equal to A^. In 
all the chapter except in subsection |3.2| the system is three-dimensional. We restrict 
to the deeply Bose-condensed regime where the non-condensed fraction is small. 
This implies that the temperature T is much lower than the critical temperature Tc 
and that the system is weakly interacting. Interactions between the cold bosons are 
characterized by the i-wave scattering length a, that we take positive for repulsive 
interactions. The microscopic details of the interaction potential are irrelevant here 
since the interaction range is much smaller than the typical de Broglie wavelength 
of the particles. The weakly interacting regime, in the considered low temperature 
regime, is then defined by [pa^Y''^ <^ 1 where p = N/V is the mean density. 

We assume that the gas is prepared in thermal equilibrium at negative times with 
some unspecified experimental procedure generally implying a coupling with the 
outer world. For clarity we consider first that the system is prepared either in the 
canonical or the microcanonical ensemble, then we apply our theory to a more gen- 
eral ensemble: a statistical mixture of microcanonical ensembles with weak rela- 
tive energy fluctuations. After the preparation phase, at positive times, the system 
is supposed to be totally isolated in its evolution. This implies that the total parti- 
cle number A^ and the total energy E are exactly conserved in time evolution. This 
assumption is realistic for ultra-cold atom experiments: the atoms are hold in con- 
servative immaterial traps and the three-body loss rates are very low in the weak 
density limit. As we shall see, this has important consequences for the temporal 
coherence of the gas. 

A first property that we discuss in this chapter is the spatial coherence of the gas. 
This is determined by the first-order coherence function 

gi(r) = (v/^(r)v>(0)) (1) 

where the bosonic field operator ^f{r) annihilates a particle in position r. The gi 
function has been measured using atomic interferometric techniques |4|. In the ther- 
modynamic limit, gi (r) tends to the condensate density po > at large distances r. 
One refers to this property as long-range order. 

A second, more subtle property, that we discuss in detail is the temporal coher- 
ence of the gas. We define the temporal coherence function of the condensate as 

(flJ(0«o(0)) (2) 

where oq is the annihilation operator in the condensate mode that is the plane wave 
with k = 0. Contrarily to the case of gi, here the operators appear in the Heisen- 
berg picture at different times. The temporal coherence function of the condensate 
is measurable (as we argue in subsection l4.1l i but it was not measured yet. The clos- 
est analog that has been measured is the relative coherence of two condensates in 
different external or internal states at equal times ||5] |6]. The coherence time of the 
condensate is simply the half width of the temporal coherence function. Remarkably 
at zero temperature it was shown that the coherence function does not decay at long 
times, it rather oscillates Q 
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(flJ(f)flo(O)) ^ («o>e'''(^=°"/'' (3) 

where («o) is the mean number of particles in the condensate and ii{T — 0) is the 
ground state chemical potential of the gas. This implies an infinite coherence time. 
At finite temperature however one expects a finite coherence time for a finite size 
system. We find that in the thermodynamic limit this coherence time diverges with a 
scaling with the system volume V that depends on the statistical ensemble in which 
the system is prepared. 

This chapter is based on our works ||8]|9] [T0|H It is organized as follows. We give 
a pedagogical presentation of the number conserving Bogoliubov theory, a central 
tool for our problem, in section |2] We apply this theory to the spatial coherence in 
section [3] The more involved issue of temporal coherence is treated in section |4] 
In subsection l4.1l we discuss how to measure (aJ(f)ao(O)) with cold atoms. General 
considerations are given in 14. 21 showing the central role of condensate phase-change 
spreading, that is then studied for different initial states of the gas. First for a single- 
mode model in 14. Bl and for the canonical ensemble 14.41 where one of the conserved 
quantities (the particle number A^ or the energy E) has fluctuations in the initial state. 
Then for the microcanonical ensemble l4.5l where none of these conserved quantities 
fluctuates. Finally in the already mentioned more general statistical ensemble within 
a unified theoretical framework in subsectionl4.6l 



2 Reminder of Bogoliubov theory 

The central result of Bogoliubov theory IITTI is that our system can be described 
as an ensemble of weakly interacting quasi-particles. The necessity to go from a 
particle to a quasi-particle picture to obtain weakly interacting objects is due to the 
presence of the condensate that provides a large bosonic enhancement of particle 
scattering processes in and out of the condensate mode. In the initial work of Bo- 
goliubov the quasi-particles are non-interacting. We will need to include the inter- 
actions among quasi-particles that give them a finite lifetime through the so-called 
Beliaev-Landau mechanism ll7l [T2ll . Here we present a powerful formulation of Bo- 
goliubov ideas introducing the phase operator 9 for the condensate mode ifTSl : in 
addition to making the theory number conserving |1T4| [TS] [T6l , 9 will play a crucial 
role for the study of temporal coherence. 



2.1 Lattice model Hamiltonian: 

Commonly a zero range delta potential Vn = gS{ri — r2) is used to model particle 
interactions with an effective coupling constant 



' Particle losses are not discussed in this chapter Their effect on temporal coherence is weak at 
relevant times as explicitly shown in 1101 for one-body losses in the canonical ensemble. 
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(4) 



(here the i-wave scattering length is a > and m is the mass of a particle). This 
however does not lead to a mathematically well defined Hamiltonian problem, even 
for two particles. As explained in ifTTI a convenient way to regularize the theory 
while keeping the simplicity of contact interactions is to discretize the coordinate 
space on a cubic lattice with lattice spacing b. This automatically introduces a cut- 
off in momentum space, since single particle wave vectors are restricted to the first 
Brillouin zone (FBZ) of the lattice [-f , f )^. Then 

Vi2=^o%^ (5) 

where now 5 is a discrete Kronecker 5. The bare coupling constant gQ is adjusted 
to reproduce the true i-wave scattering length on the lattice iflTl . 

'' = Y^b ^^^ 

where C = 2.442749 ... is a numerical constantly The Bogoliubov method is appli- 
cable when the zero energy scattering problem is treatable in the Bom regime ifTSI 
which requires here that a^bAn this limit gQ ~ g. For the lattice model to well 
describe continuous space physics the lattice spacing b should be smaller than the 
macroscopic length scales ^ and X of the gas. The healing length ^ is defined as 



h^ 



2mE, ^ 
and the thermal de Broghe wavelength as 



P8 (7) 



27tT? 
rnksT 



(8) 



Note that in the weakly interacting and degenerate limit one has ^':$> a and A ^ a. 
The system Hamiltonian in second quantized form is 



H = Y^b' 



\jf''ho\if+Y ¥''¥¥¥ 



(9) 



where ho is the one-body hamiltonian reduced here to the kinetic energy term, ho = 
— ^\, with a discrete laplacian reproducing the free wave dispersion relation E^- = 
h k^/2m when applied over a plane wave. The bosonic field operator \(f{r) obeys 
the discrete commutation relation 



• This results from the formula g^^ ~ g ' — /pB2 



(2ir)3 tp-e- 
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[V/(ri),rnr2)] = %^ (10) 



2.2 Bogoliubov expansion of the Hamiltonian 

We split the field operator into the condensate field and the non-condensed field 
i//x(r) orthogonal to the condensate wave function 0(r): 

V/(r) = 0(r)flo + V>±(r) (11) 

where aq is the annihilation operator of a particle in the condensate mode. For the 
homogeneous system that we consider, 0(r) = l/V''^. The main idea of the Bogoli- 
ubov approach is to use the fact that the non-condensed field is much smaller than 
the condensate field to expand the Hamiltonian in powers of V/x(r). This becomes 
truly operational if one succeeds in eliminating the amplitude oq of the field on the 
condensate mode. For the modulus of oq we can use the identity 

no^N-Ni_ (12) 

with N the total particle number operator, «o = floflo the condensate particle number 

A -J 4- 

operator and N± = T.t"W±W± the non-condensed particle number operator The 
elimination of the phase of do at the quantum level is more subtle, and it was not 
performed in the original work of Bogoliubov. We introduce the modulus-phase 
representation [;13i 

ao = e'^n'J^ (13) 

with the hermitian phase operator 9, conjugate to the condensate particle number: 

[ho,9]=i (14) 

It is known that the introduction of a phase operator in quantum mechanics is a 
delicate matter |fT9l . As we explain below, our formulation is not exact but it is 
extremely accurate in the present case of a highly populated condensate mode. As 
it appears from (fl4t . there is a formal analogy with the position operator x and 
the momentum operator p of a fictitious particle in one spatial dimension. For the 
fictitious particle p is the generator of spatial translations so that 

[«O,0]=/=^e'^|«o:'^> = |«o-l:0) 

where \x) represents the fictitious particle localized in position x and |no : (/>) is the 
Fock state with «o particles in the condensate mode. As a consequence the repre- 
sentation ( fT3] l of flo has the correct matrix elements in the Fock basis. The operator 
exp(/0) is a respectable unitary operator... except when the condensate mode is 
empty where one gets the meaningless result: 
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e'^\0:<j>) = \-l:<p) (15) 

This is in practice not an issue if, in the physical state of the system, the probability 
for the condensate mode to be empty is negligible. For a finite size system the prob- 
ability distribution of hq was calculated using the Bogoliubov approach and even an 
exact numerical approach 112011211 . In the thermodynamic limit we expect that the 
probability of having an empty condensate vanishes exponentially with the system 
size at r < T^.. 

In order to eliminate the condensate phase we introduce the number conserving 
operator lfT4l[T5l 

A(r)=e-'V±(r) (16) 

The success of the elimination procedure is guaranteed since the Hamiltonian con- 
serves the particle number: Injecting the splitting of the field ( fTTT i in the Hamiltonian 
and expanding, generates a series of terms in which gq appears either with Oq or with 
V/j^ (r ) . Expanding H to second order in \j/j^ and using (fT2] i we obtain the Bogoliubov 
Hamiltonian 






A^lio-fio)A+fio[U^ + l-A^^ + 2A^A 



(17) 



We have assumed that the total particle number is fixed and equal to A^ and we 
have replaced A' by A^. Still, one obtains a grand canonical ensemble for the non- 
condensed modes, with a chemical potential jUq = goP ■ The condensate indeed acts 
as a reservoir of particles for the non-condensed modes. The expression /io = goP 
is in fact the zeroth order approximation (in the non-condensed fraction) to the gas 
chemical potential. In what follows we shall take 

IM)^gp (18) 

which is consistent with the Bogoliubov theory at this order. The terms AM in 
(fTTl) represent elastic interactions between the condensate and the non-condensed 
particles. They also appear in the simple Hartree-Fock theory. The terms A^^ and 
hermitian conjugate represent inelastic interactions where two condensate particles 
collide and are both scattered into non-condensed modes with opposite momenta. 
They are absent in the Hartree-Fock theory and they play a crucial role in explaining 
the superfluidity of the gas. 



2.3 An ideal gas of quasi-particles 

To extract the physics contained in the Bogoliubov Hamiltonian one has to identify 
the eigenmodes of the system putting the quadratic Hamiltonian in a normal form. 
We present here a brief overview, a more detailed discussion was given in lfT6ll22ll . 
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In the Heisenberg picture the equations of motion of the field operators are linear, 
provided one collects A and A^ into a single unknown; 



ihdt 



A 

At 



/lo + Mo Mo \( ^ 
-IM) -(/io + Mo)/ \A^ 



-^ 



A 

At 



(19) 



The matrix S£ is not hermitian for the usual scalar product, but it is "hermi- 
tian" for a modified scalar product of signature ( 1 , — 1 ). It has moreover a symmetry 
property ensuring that its eigenvalues come in pairs ±e,f 

We now expand the field operators over the eigenvectors of ^\ 



Mr) 
At(r) 



,/k-r 



k^^O 



Uk 
yi/2 I Vk 



(kr 



^k 



Vk 
yi/2 \Uk 



(20) 



with U^ —V^ ^ I (this is the normalization condition for the modified scalar prod- 



uct). An explicit calculation gives 



Uk + Vk 



n^P/2m 



Uk-Vk \2no + ti^k^/2m 



1/4 



(21) 



The coefficients bk and b^ obey the usual bosonic commutation relations e.g. 
[^k, bl_i] = 5k k' . Injecting the modal decomposition ( l20l i in the Bogoliubov Hamilto- 
nian (flTT i one obtains a Hamiltonian of non-interacting bosons called quasi-particles: 



(22) 



The quantity Eo{N) is the Bogoliubov approximation of the ground state energy. It 
reads 




Eo{N) 



TV 



k^O 



(23) 



In the continuous space limit bj^ — > 0, the sum over k has an ultraviolet {k — > 
oo) divergence. If one replaces go by its expression ^ expanded to first order in 
ajb, go ~ g(l + Ca/b), this exactly compensates the ultraviolet divergence and one 
recovers the Lee-Huang- Yang result 



Eo{N) 



TV 



128 



{pa') 



3n1/2 



(24) 



The Bogoliubov spectrum e^, starts linearly at low k: the quasi-particles are then 
phonons. At high k one recovers the free particle spectrum shifted upwards by /Xq: 
quasi-particles in this limit are just particles. At thermal equilibrium in the canonical 
ensemble for the original system the Bogoliubov density operator is 
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a = 



•Zboi 



.e-P"Bog with j5 = l/keT 



(25) 



where Zeog is the partition function in the Bogohubov approximation. This den- 
sity operator in the canonical ensemble for particles, corresponds in fact to a grand 
canonical ensemble, with zero chemical potential, for the quasi-particles whose 
number is not conserved. 



3 Spatial coherence 

In this section we discuss the spatial coherence properties of a weakly interacting 
Bose-condensed gas, using the Bogoliubov theory. As expected one finds long range 
order in the thermodynamic limit. To complete the discussion we briefly address the 
case of a low-dimensional system where long range order is in general lost (except 
for the 2D gas at zero temperature) but where the ideas of the Bogoliubov method 
can be adapted for quasi-condensates llT7ll23l . 



3.1 Non-condensed fraction and gi function 



In a spatially homogeneous gas, the non-condensed fraction is the ratio of the non- 
condensed density (A ^A ) and the total density p . Using the modal decomposition 
(l20b and the thermal equilibrium state (l25i . one obtains in the thermodynamic limit 
in 3D: 



{N^) (AU) 



N 



(2^ 



^ul + vi 



ii^k - 1 



-V 



(26) 



This integral has no ultraviolet (k — > oo) divergence since V^ ~ 0{l/k'^). One can 
thus take the continuous space limit Z? — > and integrate over the whole Fourier 
space. The integral has no infrared (k — )• 0) divergence either, since Ul,V^ = 
0(1 /k). In order for the Bogoliubov theory to be applicable, the non-condensed 
fraction should be small. From the result (l26t one can check that this is indeed the 
case for the degenerate pX^ ^ 1 and weakly interacting (pa^)^'^ <C 1 regime. 

The first-order coherence function ([T) in the thermodynamic limit is given in the 
Bogohubov theory by 



'(r) =p- 



d^k 

(2^ 



(1 — cosk-r) 



eP^k - 1 



-V 



(27) 



where we used the exact relation {dQ\f/j_) = 0. In the large r limit, the contribution 
of the oscillating term cos k • r vanishes and ^i tends to the condensate density. This 
implies that spatial coherence extends over the whole system size. 
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3.2 In low dimensions 



In a straightforward generalization of (126) to low dimensions, the non-condensed 
fraction is infrared divergent in 2D for T > Q, and in ID for all T: there is no Bose- 
Einstein condensate in the thermodynamic limit in agreement with the Mermin- 
Wagner-Hohenberg theorem ll24] |251 . Nevertheless, in the weakly interacting and 
degenerate regime there are weak density fluctuations and weak phase gradients. 
This is the so called quasi-condensate regime Il23ll26l . The main ideas of the Bo- 
goliubov approach can still be applied after the introduction of a modulus-phase 
representation of the field operator \jf in each lattice site 1271 : 



V/(r)=e"'Wyp^ (28) 

where p (r)fe^ and 9 (r) are conjugate variables similarly to (Il4l and d is the spa- 
tial dimension. As we discussed in subsection 12.21 and in [17), the modulus-phase 
representation of the annihilation operator in a given field mode is accurate if this 
mode has a negligible probability to be empty. This in particular requires that the 
mean number of particles per lattice site is large, pb'' ;» 1. In the weakly interacting 
p(^'' 3> 1 and degenerate pA'' ^ 1 regime, one can adjust b to satisfy this condition 
while keeping b^£,,X so as to well reproduce the continuous space physics. In this 
regime one also finds that the probability distribution of the number of particles on 
a given lattice site is strongly peaked around the mean value pb'' ^ I, with a width 
much smaller than the mean value, which legitimates the representation ( l28T l. 

If one blindly applies the plain Bogoliubov result (|27] | in the absence of a con- 
densateQ one finds that the first-order coherence function g ^ °° (r ) — > — oo at infinity, 
logarithmically with r in 2D {T > 0) and in ID {T = 0), and even linearly in r in 
lDatT> 0. One may believe at this stage that ^j °°(r) is simply meaningless in 
those cases. The extension of the Bogoliubov theory to quasi-condensates however 
produces the remarkable result ll27l : 



^?*^(r)=Pexp 



B02/ \ 

^1 ir) 



-1 



(29) 



The quasi-condensate first-order coherence function g j (r) tends to zero for r — s> oo 
as a power law in 2D (T > 0) and in ID {T = 0), and exponentially for T > in ID, 
as expected l23\. The gas has then a finite coherence length l^ (e.g. the half-width of 
gi) much larger than ^ or A in the weakly interacting and degenerate regime. Over 
distances r <C k, phase fluctuations are small, and the system gives the illusion of 
being a condensate: one can linearize the exponential in Eg. ([29b, to obtain g^ (r) ~ 
gj°^(r). The phase and density fluctuation properties of the quasi-condensates at 



^ One may wonder in 2D about the value of ^o = gop, since go logarithmically depends on the 
lattice spacing b |27|. and dimensionality reasons prevent from forming a coupling constant g 
(such that gp is an energy) from the quantities h, m and a, where a is now the 2D scattering length, 
given in I26II28I . According to 1271 one simply has to take for jlo the gas chemical potential ll{T). 
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nonzero temperature have been studied experimentally with cold atoms in ID 
and in 2D 1321 [33l and confirm the theoretical picture. 



4 Temporal coherence 

In this section we discuss the temporal coherence properties of a finite size Bose- 
condensedgas, defined by the coherence function (a,j(f )flo(0)) already introduced in 
equation (|2]i. Although, strictly speaking, this coherence function was not measured 
yet with cold atoms, we argue in section WA\ that it is in principle measurable. In 
subsection |4.2| we show that the condensate coherence function ^ can be related to 
the condensate phase-change during the time interval t. The loss of temporal coher- 
ence is thus due to the spreading in time of this phase-change, which is the quantity 
that we actually calculate. Whenever one of the conserved quantities (total particle 
number A^ or total energy E) fluctuates in the initial state from one realization to 
the other, the phase-change spreads ballistically. Once the effect of fluctuations of 
A^ is understood (subsection I4.3l l, the more involved effect of energy fluctuations 
for fixed A^ can be understood by analogy. The resulting guess for the phase-change 
spreading can be justified within the quantum ergodic theory (subsection l4.4l ). The 
only case in which pure phase diffusion is found is when the conserved quantities A^ 
and E are fixed, that is in the microcanonical ensemble (subsection l4.5l ). For fixed A^ 
and a general statistical ensemble for energy fluctuations, we finally give in subsec- 
tion |4J6] the expression for the variance of the phase-change in the long time limit, 
that includes both a ballistic term and a diffusive term. 



4.1 How to measure the temporal coherence function 

We give here an idea of how to measure the condensate temporal coherence function 
(flj(f )ao(0)) in a cold atom experiment | lOl. The scheme uses two long-lived atomic 
internal states \a) and \b) and it is a Ramsey experiment as in Q, with the notable 
difference that the pulses are arbitrarily weak instead of being k/2 pulses. 

The Bose-condensed gas is prepared in equilibrium in the internal state \a) and 
the state \b) is initially empty. At time zero one applies a very weak electromagnetic 
pulse, of negligible duration, coherently coupling the two internal states. After the 
pulse, the system evolves during a time t in presence of interactions only among 
atoms in \a): we assume no interactions between a and h component^ and neg- 
ligible interactions within the b component due to the very weak density in that 
component. At time t one applies a second pulse of the same amplitude, and one 
measures the particle number in state \b) in the plane wave k = 0. 



* This can be realized experimentally either using a Feshbach resonance 1341 or spatially separating 
the two components 1351 . 
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The scheme can be formalized as follows. The first pulse, at f = 0, coherently 
mixes the two bosonic fields !//„ and xj/j, with a real amplitude rj so that 



Va(r,0+) = ^l - ri'- Wair, 0-) + ri Whir, 0-) (30) 

Wb{r,0+) = v/l-7?2v/i(r,0-)-77V/„(r,0-) (31) 

In between time 0+ and time t^ the two fields evolve independently. Field \J/„ 
evolves in presence of kinetic and interaction terms as in (|9]l. Field 1//;, evolves with 
kinetic and internal energy terms so that its amplitude on the k = mode obeys 

bo{t-)=e'^'bo{0+) (32) 

where S is the detuning between the electromagnetic field and the a — b atomic 
transition (the calculation is peiformed in the rotating frame). The second pulse 
at time t mixes again the two fields with the same mixing amplitudes as in (l30l l. 
([3TI 1. After the second pulse one measures Nho{t) ~ {{p]p^'){t^')) . Using the mixing 
relations and ( |32] | one expresses /?o(f^) as a function of/7o(0^),flo(0^) and ao(f )■ 
Since the initial state for component b is the vacuum, the contribution of ^o(O^) 
vanishes and one obtains the exact relation: 

A^i0(0 = T7'{(l-77^)((4«0)(0")) + ((4«0)(f"))pulse 

+ v/T^^[e'*(4(f-)flo(0-))puLse + c.c.]} (33) 

that we expand for vanishing r\: 

A^io(0=2T?'{(«o)+Re [e'*(fl5(r)flo(0))]} +0(774) (34) 

In particular, the subscript (. . .)puise on the expectation values, indicating that they 
are taken for a system having experienced the first pulse, was removeq5 The desired 
correlation function (aJ(f)ao(O)) can be extracted from the contrast of the fringes 
obtained by varying the electromagnetic field frequency. The signal A^foo(0 itself is 
small (it is proportional to 77^) but the contrast of the fringes is independent of 77 in 
the small 77 Umit, and it starts at unity at r = 0. 



^ The expectation values (. . .)puise differ from the original ones (. . .) in the absence of pulse by 
0(r\-\. To first order in r\, the perturbation of !//„ due to the pulse is linear in i/i'/,(0^) and has a 
zero contribution to the expectation values since component b is initially in vacuum. 
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4.2 General considerations about (aQ(?)ao(0)) 

4.2.1 Phase-change spreading 

Here we go through a sequence of transformations that relates the temporal co- 
herence function (aQ(f)ao(O)) to the variance of the condensate phase-change 
9{t) — 0(0). We use the modulus-phase representation (fTjt of the annihilation op- 
erator flo- Since the non-condensed fraction is very small, we simply neglect the 
fluctuations of the modulus of aq ie- we replace «o with its mean value in equation 
([H. We then obtain^ 

(4(f)flo(0))^(«o)(e-'[^"'-^'''") (35) 

If the phase-change (r) — (0) has a Gaussian distribution, which may be checked 
a posteriori, the application of Wick's theorem yields 



7j(f)flo(0)) ~ //5„\.-'(e(0-e(0)>.-Var[e(r)-e(0)]/2 



(no)e ^ " ' "e 



(36) 



This remarkable formula quantitatively relates the loss of temporal coherence in an 
isolated Bose-condensed gas to the spreading of the condensate phase-change. 

The operational way to determine the condensate phase-change spreading is to 
work with the phase derivative: contrarily to 0, is a single-valued hsrcnitmn oper- 
ator that has a simple expression within the Bogoliubov approach. The correlation 
function of the phase derivative 

C(0 = (^(0^(0))-(^)' (37) 

gives access to the variance of the phase-change by simple integration: 



Var [0(f)- 0(0)] =2r / dxCR{x)-2 [ dxxCR{x) 



(38) 



where Cr is the real part of C. One obtains a single integral (rather than a double 

integral) using the fact that the real part of {Q{t\)Q{t2)) is a function of |fi —to] 
only, for a system at equilibrium. The long-time behavior of Cr determines how the 
phase-change spreads at long times as summarized in Fig.[T] 

At finite temperature, one might expect that 0(f) decorrelates from 0(0) at long 
times so that Cr ^ Q and the phase-change spreading is diffusive. As we will see, 
this is however not the case, except if the system is prepared in the microcanonical 
ensemble. This is a consequence of energy conservation between times and f in our 



* Here we have neglected the non-commutation of 9(f) and 9(0). From the Baker-Campbell- 
Hausdorff formula, and to zeroth order in the non-condensed fraction, see equation | |45K the cor- 
rection is a factor e^ 21'' ('V)+''('V ') which is in'elevant for our discussion. 
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isolated system. This point was overlooked in the early studies of ||36l[37l[38ll where 
the non-condensed modes were treated as a Markovian reservoir and phase diffu- 
sion was predicted. A subsequent study ||39l based on a many-body Hamiltonian 
approach showed that phase-change spreading is ballistic for a system prepared in 
the canonical ensemble. The coefficient of f ^ in ||39l was however calculated within 
the pure Bogoliubov approximation, neglecting the interactions between the Bo- 
goliubov quasi-particles, which is illegitimate in the long time Umit as we shall see. 



diffusive regime 



Cr{x) ^=^J,{\/X) 



Var[e(0-9(0)] ~ 2Df 



Cr(^) 



lVar[9(t)-e(0)]~2ti;dxq(T) 



ballistic regime | 


lim^^+oo Ce 


(t)=A 


7^0 


Var[e(f)- 


6(0)] ~ 


At' 




X 

Fig. 1 Different regimes of the condensate phase-change spreading at long times. Cr is the real 
part of the congelation function C defined in {37j. 



4.2.2 Key ingredients of the theory 



In order to correctly determine the phase-change spreading in the long time limit, we 
shall use two key ingredients in our theoretical treatment: an accurate expression of 
the phase derivative and the inclusion of the interactions among Bogoliubov quasi- 
particles, to which we add the constraint of strict energy conservation during the 
system evolution. 

Time derivative of condensate phase operator. The commutator of with H given 
by ^ is calculated exactly using 



[^:V^(r)] = -ao 



'•0(r) 
2«o 



(39) 



and its hermitian conjugate, with the condensate wave function 0(r) = 1/V'/^. The 
exact result is given in equation (67) of |8j|. Expanding up to second order in the 
non-condensed field A and using the modal decomposition (|20|) . one obtains for 
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fixed Af:|Zl 



l[e,H].-i,(r^o)-fL(^. + nf% 



(40) 



We have introduced the zero-temperature chemical potential /i(r = 0) = -^Eq[N), 
where Eq{N) is given in ( |23] ). and the quasi-particle number operators 

«k = blh (41) 

The expression (|40] | of the phase derivative differs from the one heuristically intro- 
duced in ll37l[38l : 9 is not simply equal to —giiQ/hV. 

Interactions between quasi-particles: Pushing one step further the Bogoliubov ex- 
pansion of section |2] that is including terms up to third order in the non-condensed 
field, one obtains 

Hc^Hjiog + H3 (42) 

where ^Bog is the Bogoliubov Hamiltonian ( l22b and 



^3=^op'/'L^'a + (A+A^)A (43) 

r 

The Hamiltonian Hj, is cubic in the field A and it corresponds to interactions be- 
tween quasi-particles. While Hboi> is integrable (all the i^ are conserved quantities), 
the Hamiltonian //eog +^3 is not integrable, which plays a central role in condensate 
dephasing. By replacing A with its modal decomposition ( l20l i in i/3, two types of 
resonant processes appear, that do not conserve the total number of quasi-particles: 
the b^b^b Beliaev process and the b'bb Landau process. In the Beliaev process 
one quasi-particle decays into two quasi-particles, while in the Landau process two 
quasi-particles merge into another quasi-particle. The processes involving b^b^b^ 
and bbb are non-resonant (they do not conserve the Bogoliubov energy) and they 
cannot induce real transitions at the present order 



4.3 If N fluctuates 

In this subsection we allow fluctuations of the total number of particles and we 
investigate their effect on temporal coherence. The effect is akeady present in the 
case of a pure condensate, so that we restrict to a one-mode model in this subsection: 
identifying the condensate particle number «o with the total particle number N, we 
obtain the model Hamiltonian 



' We have neglected oscillating tenns in bb and wb'': after time integration of 9 they give a 
negligible contribution to 6{t) — 9(0). 
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'^one mode — or/ \^) 

The condensate phase derivative is 

B{t) = ^[6, Hon, mode] = -^{N)/n (45) 

in 

where the chemical potential for the system with A^ particles is simply ji (N) = gN/V 
for the one-mode model. Since A' is a constant of motion, temporal integration is 
straightforward: 

e{t)-e{o) = -ii{N)t/h (46) 

If A^ is fixed there is no phase-change spreading. If the initial state is prepared with 
fluctuations in A^ then the phase-change spreads ballistically 114011411 : 



Var[0(O-0(O)] = (fAzr(^^j VarA^ (47) 

Correspondingly the temporal coherence function {0^(1)00) decays as a Gaussian 
in tun g^ll42l |43l . A similar phenomenon was observed experimentally B4l|45]|46l 
not for the temporal correlation of a single condensate but for equal-time coherence 
{aQ{t)bo{t)) between two condensates prepared in different modes or internal states 
with a well defined relative phase and fluctuations in the relative particle number 



4.4 N fixed, E fluctuates: Canonical ensemble 

We assume in this subsection that the gas is prepared in equilibrium at finite tem- 
perature T in the canonical ensemble with A^ particles. We first treat this case by 
analogy with the previous subsection, and then we expose a systematic derivation 
of the result based on quantum ergodicity. 



4.4.1 Using an analogy with the case of fluctuating A' 

Similarly to A' in the previous subsection, here H is a conserved quantity that fluc- 
tuates in the initial state. Indeed the canonical ensemble is a statistical mixture of 
energy eigenstates with different eigenenergies. By analogy with ( |46] | we expect that 

9it)~9(p)^-^i^,iH)t/h (48) 

where /imc {E) is the chemical potential of the microcanonical ensemble of energy E. 
As relative energy fluctuations are vanishingly small for a large system, we can lin- 



The phase revivals at macroscopic times multiples of iTttiV /g I42II43I are absent here due to the 
Gaussian hypothesis used to obtain <36b . 
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earize Hmc{E) around the mean energy E to obtain a ballistic phase-change spread- 
ing 



Var[0(r)-0(O)]-(r/ri)2 



dE ^ ' 



Var// (49) 



The coefficient of t^ is proportional to the variance of the energy in the initial state 
and scales as the inverse of the system volume in the thermodynamic limit. For con- 
venience, one can reexpress this coefficient in terms of canonical ensemble quanti- 
ties using UraclEiT)] = niT) (for a large system) so that j^HmdE) = jj,fi/-^E, 
where /i(r) and E{T) are the chemical potential and mean energy in the canonical 
ensemble at temperature T. An explicit expression of the coefficient of f^ is given in 
Eq. (73) of |8| using Bogoliubov theory to evaluate the partition function, E{T) and 
li{T). The obtained formula for /i(r) also gives the intuitive and interesting side 
result 

{e)^-pi{j)in (50) 



4.4.2 From quantum ergodic theory 

In the previous analogy leading to (|49] l there is a strong implicit hypothesis. The 
fact that the phase-change is a function of the Hamiltonian only, see Eq. ( |48] l, is in 
general true only for an ergodic system in the long time limit. For example if the 
Hamiltonian was truly equal to //eog- ^(0 ~ ^(0) would depend on the set of all 
occupation number operators i\ and Eqs. ( 148149b would not apply. 

We now derive Eq. ( |49] l using quantum ergodic theory. To this end we calculate 
the asymptotic value of the correlation function C{t). To eliminate oscillations of 
C{t) we evaluate its time average. By inserting a closure relation over exact eigen- 
states If;^) with eigenenergies Ex of the interacting many-body system, we obtain 

\ fdzCiz) -> i:p;t|(f^|^|f^)|^- (E/^;i(f;il^|f;i) 1 (51) 

where /?;t is the probability to find the system in the eigenstate | f^ ) . In the canonical 
ensemble p^ = e\p{—l5Ex)/Z. In ( BTI ) we have assumed that there are no degenera- 
cies consistently with the non-integrability of the systerrQ For a classical system, 
ergodicity implies that the time average over a trajectory of energy E coincides with 
the microcanonical average at that energy. The extension of this concept to a quan- 
tum system is the so-called eigenstate thermalization hypothesis 1471 l48l |49l : the 
mean value of a few-body observable (9 in a single eigenstate |f^) is very close to 
the microcanonical average at the same energy: 



' For a large system the level-spacing 5E vanishes exponentially with the system size, and one 
may fear that an exponentially long time / > h/SE is needed to reach the limit J5U . However, the 

corresponding off-diagonal matrix elements of 9 also vanish exponentially with the system size in 
the eigenstate thermalization hypothesis 1471 . 
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(f;i|0|fi)~4nc(£ = £;i) (52) 

We apply this hypothesis to the operator O ^ 0. The last step is to realize that 

within the Bogoliubov theory, the microcanonical average of 6 is proportional to 
the microcanonical chemical potential I "l 

0n,e(£) = -Mmc(£)//J (53) 

One then obtains 



f2 



Var[0(r)-0(O)] - -.Wav^^^iH) 
f-1.00 jj-^ 



(54) 



Linearizing /Zmc (H) in ( l54l i for small relative energy fluctuations around E one re- 
covers (|49] |. 



4.4.3 Physical implications 

A consequence of (l49T l is that, for a system prepared in the canonical ensemble, 
the correlation function C(t) of 9 does not tend to zero when T -^ +°°. The same 
conclusion is reached for the correlation function of 1%, whose long time limit can 
be calculated with the quantum ergodic theory |i8J. This qualitatively contradicts 
ll36l l37l [38l ■ It only qualitatively agrees with ||39l since the system Hamiltonian H 
in ||39l was eventually replaced by the integrable Hamiltonian Hso„. 

In ll36l [37I [38JI the non-condensed modes were treated as a Markovian reservoir 
This approximation is excellent to calculate temporal correlation functions of "mi- 
croscopic" observables such as the quasiparticle numbers. For example, this gives 
fork,k'7^0||8l: 

(nk(0«k'(0)) - («k)(«k') ^=°' 4,k'(«k)(l + («k))e-^''' (55) 

where the damping rate fk is due to the Beliaev-Landau processes. However quan- 
tum ergodic theory shows that the exact long time limit of this correlation function 
is nonzero (even for k ^ k') but rather a quantity of order 1 /N. In the double sum 
over k and k' that appears in C(t), this introduces a macroscopic correction of order 
A^ missed by the Markovian approximation. 

We illustrate this discussion in Fig|2]with a classical field model lID. The ex- 
act numerical result (black squares linked by a solid line) confirms the ergodic re- 
sult (dash-dot-dotted blue curve). The flat red dashed line is the Bogoliubov theory 
where the «k are constants of motion. It is close to the numerical result only at short 
times. The dash-dotted violet curve that tends rapidly to zero is a Markovian model 
based on ( fSSl l. 



"^ See reference [45] of (8]. In fact for a large system it is sufficient to prove the equality in the 
canonical ensemble of mean energy E, as already given by Eq. l lSOt . 
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Fig. 2 For a gas prepared 
in the canonical ensemble, 
correlation function of 9 for 
the classical field model. The 
equation of motion is the non- 
linear Schrodinger equation. 
This coiTesponds to Fig. 6 in 
tSl . V is the volume. See text 
for the meaning of the various 
curves and symbols. 
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4.5 N fixed, E fixed: Microcanonical ensemble 



In this section we assume that the gas is prepared in the microcanonical ensemble of 
energy E. According to ( |54l i the coefficient of the ballistic spreading of the phase- 
change is zero. It was found in ifTOl that C(t) = 0(1/t^) at long times, so that the 
phase-change spreads diffusively, with a diffusion coefficient defined by 



Var[0(f)-0(O)]-2Df with 



£)= / dxCRix) 
/o 



(56) 



To determine D we thus need the whole time dependence of C(t). From (|40l i. C{x) 
can be deduced from all the correlation functions (nk('r)«k'(0)) of the quasi-particle 
number operators. Within the Bogoliubov approximation for the initial equilibrium 
state, the gas is prepared in a statistical mixture of Fock states !{«„}) of quasi- 
particles where, in any given Bogoliubov mode of wave vector q, there are exactly 
Hq quasi-particles («q is an integer). One can then calculate the correlation functions 
for an initial Fock state !{«„}) and average over the microcanonical probability 
distribution for the {«„}. 



For a given initial Fock state, one then simply needs 

„k(T)^(K}|«u(T)|K}) 



(57) 



In the thermodynamic limit, the evolution of such mean numbers of quasi-particles 
are given by quantum kinetic equations including the Beliaev-Landau processes due 
toi/3 1501: 

nk^--j—2 M^?[nknq-nk+q(l+«q + nk)] [f^lk) 5(e;t + £,, - e|k+q|) 
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with the Behaev-Landau couphng ampUtudes: 



<i' 



^ UciUkUk> + V,ykVk' + {Uci + V,) {VkUy + UkVk> 



(59) 



The first Hne in (|58] | describes Landau processes and the second Hne describes BeH- 
aev processes. In practice we Hnearize the kinetic equation (l58l l around the equihb- 
rium solution «k LJ and we solve the resulting linear system numerically. We refer 
to ifTOl for technical details. 

The phase diffusion coefficient is shown in Fig|3]as a function of the temperature 
T such that the mean canonical energy E{T) is equal to the microcanonical energy 
E. Remarkably, when D and T are properly rescaled (as in the figure), the curve 
is universal. In particular this shows that D vanishes as the inverse of the system 
volume in the thermodynamic limit. Interestingly, at low temperature, D vanishes 
with the same power-law T^ as the normal fraction of the gas: 



tiDV 



, 03036 ('^X 



(60) 



We performed classical field simulations in the microcanonical ensemble fpl. 
As expected we found that the phase-change has a diffusive behavior: its variance 
increases linearly in time at long times (not shown) and the phase-change probability 
distribution is well adjusted by a Gaussian as we show in the left panel of Fig. H) In 
the right panel Fig. |4]we show that the diffusion coefficient is well reproduced by a 
classical field version of the kinetic theory. 
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' ' For an infinite system, the stationary solution of | |58) is ensemble independent and corresponds 
to the Hose formula Tx\^(E) = l/(expj3et- — 1), where j8 is adjusted to give the mean energy E. 
Finite size effects on the h^, that can be calculated from Eq. (61) of fS), are here not relevant. 
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4.6 A general statistical ensemble 

We now consider a generalized ensemble at fixed A^ that includes both the micro- 
canonical and the canonical ensembles as particular cases. This is a statistical mix- 
ture of microcanonical ensembles with a probability distribution P(£') of the system 
energy E that depends on the particular experimental procedure to prepare the ini- 
tial state of the gas. Remarkably the approach of the previous subsection based on 
kinetic equations can be extended to this case. 



4.6.1 General result for the phase-change spreading 

Provided that the relative energy fluctuations vanish in the thermodynamic limit, we 
find the long time limit |ITOl 




(61) 



For the coefficient A of the ballistic t^ term we recover the form of the quantum 
ergodic result (|49] l. This is not surprising as the reasoning of subsection l4.4l does not 
rely on the fact that the system is prepared in the canonical ensemble. On the other 
hand the value of the coefficient does depend on the statistical ensemble through the 
mean energy E and the variance of the energy. A physical derivation of this result 
within kinetic theory is given in the next subsection. 

A remarkable result is that, in the general ensemble, the phase derivative correla- 
tion function C(t) is the sum of its long time limit A and of the correlation function 
Cmc(f ) in the microcanonical ensemble of energy E: 




9 12 15 IS 21 24 

Fig. 4 Classical field simulations in the microcanonical ensemble. Left panel (taken from (9|): 
Probability distribution P{<p) of the condensate phase-change fluctuations (j) = 9(f) — 6(0) — 
{6(/) — 9(0)) at a large time t. The dashed line is the expected Gaussian. Right panel (taken from 
1101 ): Diffusion coefficient as a function of the temperature, extracted from the numerics (bullets 
with error bars) and calculated by the classical field version of the kinetic equations l l58t (crosses 
linked by segments). 
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C(T)=A + Cn,c(T) (62) 

As a consequence the diffusion coefficient D of Eq. dHTT i is the same as the one 
for the microcanonical ensemble of energy E. The same conclusion holds for the 
constant time offset foffLj: 



D= / dTCR.mciT) (63) 

Jo 

/o°°t/TTCfl,mc(T) 

foff = r°°^^/- 77V ^^^> 

JQ flTQ,mc(T) 

where C/;,mc is the real part of Cmc- The physical origin of the time offset foff is 
apparent in Eq.(l64li: it is due to the finite width of the phase derivative correlation 
function. As CR_mc{'^) is found to be positive, foff can be simply interpreted as the 
coiTelation time of the phase derivative in the microcanonical ensemble. The formal 
expressions for D and foff, in terms of the matrix of the linearized kinetic equations, 
are given in ifTOI . 

These results are made more concrete by Fig. |5] for a quantum system in the 
thermodynamic limit, we show the microcanonical correlation function Cmc(0 as 
a function of time, and the variance of the phase-change either in the canonical 
ensemble of temperature kgT = lOpg or in the microcanonical ensemble with the 
same mean energy. This reveals in particular that the asymptotic expression dHTT l 
becomes rapidly accurate. 



4.6.2 Recovering the ballistic spreading from kinetic theory 

Due to energy conservation, the linearized kinetic equations have a zero-frequency 
undamped mode. We will show that, in presence of energy fluctuations in the initial 
state, the amplitude over this mode is nonzero, so that the phase derivative correla- 
tion function C(t) does not tend to zero at long times and the phase-change variance 
shows a t^ term as in Eq. (1611 1. The derivation presented here was significantly sim- 
plified with respect to the original one of ifTOl . 
We introduce the notation 

«k(£')=nkmc(£') (65) 

for the average number of quasi-particles in mode k in the microcanonical ensemble 
of energy E. The kinetic equations ( fSSl l. linearized around the stationary solution 
{«q(£)}, can be put in the form 

i(T)=Mx(T) (66) 



^^ This is true to leading order in the system size since our linearized kinetic approach cannot 
access the subleading tenns. 
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Fig. 5 For a quantum system in the thennodynamic limit, the microcanonical phase derivative cor- 
relation function Cmc(0 (red solid line, light veilical axis) and the variance of the phase-change 
(black lines, left vertical axis) are shown as functions of time. For the variance, the upper (lower) 
solid line is for the canonical (microcanonical) ensemble, and the dashed lines are the correspond- 
ing asymptotic fonns of Eq. | |6U . kgT = lOpg, V is the system volume, g is the effective coupling 
constant (|4) and ^ is the healing length Q. This is Fig. 3 of |10| . In atomic condensates £, is in the 
;Um range and the time unit of the figure is in the ms range. 



where we have collected all the unknowns nk(T) — nk{E) in a single vector x(t) and 
M is a matrix. The existence of a zero frequency mode can be understood in two 
different ways that we explain. 

First reasoning: Consider an energy E close to E. In the same way as {nk(£)}, the 
set of occupation numbers {n]i{E)} constitutes a stationary solution of the full ki- 
netic equations ( fSSl l. Since the solutions are close, their difference {«k(£') — «k(£')} 
obeys the linear system (|66] | so that the vector eo of components 



eo,k = ^"kl^) 



(67) 



is a zero-frequency eigenmode of M. 

Second reasoning: The Bogoliubov energy £k^o £*:«k ('^) is conserved by the kinetic 
equations. An a consequence e ■ x(t) is a constant (the vector e has components 
e^,) and its time derivative is zero. This holds for all initial values of x , and thus 
implies that e is a left eigenvector of M with zero eigenvalue. A basic theorem 
of linear algebra then implies the existence of a right eigenvector of M with zero 
eigenvalue. Actually we already found it: it is eo of components ( |67] |. Such left 
and right eigenvectors are called adjoint vectors. For our normalization choice, their 
scalar product e • eo = jgE = 1 as it should be. 

We now go back to the correlation function C(t). We introduce the (zero-mean) 
fluctuation operators 

5«k = nk-nk(£) (68) 
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where we have neglected the difference between («k) and nk{E) in the large system 
size limit. The correlation function C(t) is then obtained as 

C{t)^A-x{t) with ;ck(T) = -(5«k(T)^(0)) (69) 

where we have collected in a vector A, the coefficients in given by Eq. ( l40b : 

Ak = ^iUk + Vkf (70) 

Following the reasoning of subsection l4.5l on finds that x{t) obeys Eq. ( |66] ). Splitting 
x{t) = yeo +X (t) we have in the long time hmit that X (t) -^Q due to the Beliaev- 
Landau damping processe whereas 7= e •.y(0) is a constant. At long times one then 
has 

C(t) ^ [e-mKA-eo) (71) 

Taking the microcanonical average of (|40] | and using (|53) on obtains the Bogoliubov 
expression for the microcanonical chemical potential: 

Mmc(£) =^(7^ = 0)^+^ ftAk«k(£) (72) 

kT^O 

Using the expression of eo this leads to A • eo = -j^p^mc{E)/h. We now evaluate the 
expectation value (. . .} appearing in e ■x{Q) in two steps. We first take the expec- 
tation value in the microcanonical ensemble of energy E: one can then replace the 
operator £k £Ar5''k(0) with E — E, since the total Bogoliubov energy is fixed to E. 
One is left with a microcanonical average of 0(0) at energy E, an average already 
given by Eq. ( |53] ). and that one can expand around E to first order in E ~ E. The 
last step is to average over E with the probability distribution P{E) defining the 
ensemble, to obtain 

e.x(0)== Var(£)^(£) (73) 

Collecting all the results, we exactly recover the coefficient of f^ in Eq. dMT l. 

After this last reasoning, it becomes apparent that, contrarily to the zero-frequency 
component jeQ, the contribution of the damped component X{t) of x(t) can be 
treated to zeroth order in the energy fluctuations: one can directly take E ^ E with- 
out getting a vanishing contribution to C(t) and to Eq. ( |6TI ). This explains why both 
the diffusion coefficient D and the time offset foff, that purely originate from X{t), 
are essentially ensemble independent. 
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